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Abstract 

The CLEO experimental data on the tt'j transition are analyzed to next-to-leading order accuracy 
in QCD perturbation theory using light-cone QCD sum rules. By processing these data along the 
lines proposed by Schmedding and Yakovlev, and recently revised by us, we obtain new constraints 
for the Gegenbauer coefficients 02 and £14, as well as for the inverse moment (x -1 )^ of the pion 
distribution amplitude (DA). The former determine the pion DA at low momentum scale, the 
latter is crucial in calculating pion form factors. From the results of our analysis we conclude 
that the data confirm the end-point suppressed shape of the pion DA we previously obtained with 
QCD sum rules and nonlocal condensates, while the exclusion of both the asymptotic and the 
Chernyak-Zhitnitsky DAs is reinforced at the 3a- and 4<r-level, respectively. The reliability of the 
main results of our updated CLEO data analysis is demonstrated. Our pion DA is checked against 
the di-jets data from the E791 experiment, providing credible evidence for our results far more 
broadly. 
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1. Introduction 



The recent high-precision CLEO results [lflfor the 7T7 transition form factor gave rise to 
dedicated theoretical investigations jE IE 0HE IE IE IE 13 • These experimental data are of 
particular importance because they can provide crucial quantitative information on nonper- 
turbative parameters of the pion DA and — as we pointed out in [9] — on the QCD vacuum 
nonlocality parameter A 2 , which specifies the average virtuality of the vacuum quarks. In 
the absence of a direct solution of the nonperturbative sector of QCD, we are actually forced 
to extract related information from the data, relying upon a theoretical analysis as complete 
and as accurate as currently possible. 

It was shown by Khodjamirian |E] that the light-cone QCD sum- rule (LCSR) method 
provides the possibility to avoid the problem of the photon long-distance interaction (i.e., 
when a photon goes on mass shell) in the 7*(Q 2 )7(<Z 2 ) — ► tt° form factor by performing 
all calculations for sufficiently large q 2 and analytically continuing the results to the limit 
q 2 = 0. Schmedding and Yakovlev (SY) jE| applied these LCSRs to the next-to-leading order 
(NLO) of QCD perturbation theory. More recently |9j, we have taken up this sort of data 
processi ng in an attempt to (i) account for a correct Efremov-Radyushkin-Brodsky-Lepage 
(ERBL) [lQj evolution of the pion DA to every measured momentum scale, (ii) estimate more 
precisely the contribution of the (next) twist-4 term, and (iii) improve the error estimates 
in determining the la- and 2<r-error contours. 

The main outcome of these theoretical analyses can be summarized as follows: 



the asymptotic pion DA [10| and the Chernyak-Zhitnitsky (CZ) [llj model are both 
outside the 2<r-error region 

the extracted parameters and (i.e., the Gegenbauer coefficients of the pion DA) 
are rather sensitive to the strong radiative corrections and to the size of the twist-4 
contribution 

the CLEO data allow us to estimate the correlation scale in the QCD vacuum, A 2 , to 
be < 0.4 GeV 2 . 



The present note gives a summary of our lengthy analysis |9j extending it a few steps 
further, notably, by obtaining from the CLEO data also a direct estimate for the inverse 
moment of the pion DA that plays a crucial role in electromagnetic or transition form factors 
of the pion and by verifying the reliability of the main results of the CLEO data analysis 
quantitatively. Moreover, we refine our error analysis by taking into account the variation of 
the twist-4 contribution and treat the threshold effects in the strong running coupling more 
accurately. The predictive power of our updated analysis lies in the fact that the value of 
the inverse moment obtained from an independent QCD sum rule is compatible with that 
extracted from the CLEO data, referring in both cases to the same low-momentum scale 
of order of 1 GeV. As a result, the pion DA obtained before 0] from QCD sum rules with 
nonlocal condensates turns out to be within the la-error region, while the asymptotic and 
the CZ pion DAs are excluded at the 3a- and 4<r-level, respectively. Our predictions for the 
pion DA are found to be in agreement with the Fermilab E791 data 2 • 
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2. Light cone sum rules 



Below, we sketch the improved NLO procedure for the data processing, developed in 
0. Let us recall that this procedure is based upon LCSRs for the transition form factor 
F 7 * 77r (Q 2 , q 2 « 0) (HEf. Accordingly, the main LCSR expression for the form factor 



SQ oo 



F?dW 2 )= 1 - j^ P {Q\s-y?)eH-s)/^ + l - |^p(Q 2 ,*;p 2 ) (2.1) 



p 

s 



follows from a dispersion relation. The corresponding spectral density p(Q 1 , s; p ) = 



I m 



F^(Q 2 ,q 2 = -s;p 2 ) 



is calculated by virtue of the factorization theorem for the 

form factor at Euclidean photon virtualities q\ = —Q 2 < 0, q\ — —q 2 < [HI EH, 
with M 2 « 0.7 GeV 2 being the Borel parameter, whereas m p is the p-meson mass, and 
so = 1.5 GeV 2 denotes the effective threshold in the p-meson channel. The factorization 
scale p 2 was fixed by SY at p 2 = p§ Y = 5.76 GeV 2 . Moreover, Fq ( j D 7T (Q 2 , q 2 ; p 2 ) contains a 
twist-4 contribution, which is proportional to the coupling £ 2 (p 2 ), defined bv pi 13 



(n(p)\g s dG aill a u\0} = ib'UVy, , (2.2) 

where G ap = {l/2)e aixp(T G^ and G pa = G a pa \ a /2. 

This contribution for the asymptotic twist-4 DAs of the pion as well as explicit expressions 
for the spectral density p(Q 2 ,s;p 2 ) in LO have been obtained in [Hf to which we refer for 
details. The spectral density of the twist-2 part in NLO has been calculated in [f| — see Eqs. 
(18) and (19) there. All needed expressions for the evaluation of Eq. (J2.1)) are collected in 
the Appendix E of [9], cf. Eqs. (E.1)-(E.3). 

We set p 2 = Q 2 in Fqqy) 7 \Q 2 i Q 2 '-! an d use the complete 2-loop expression for the form 
factor, absorbing the logarithms into the coupling constant and the pion DA evolution at 

the NLO level P so that a s (p 2 ) a s (Q 2 ) (RG denotes the renormalization group) and 

^(xsp 2 ) ip n {x;Q 2 ) = U{fi 2 -* Q 2 )f n {x;fi 2 ). 

Then, we use the spectral density p(Q 2 , s), derived in p at p 2 = pg Y , in Eq. (|2.1j) to obtain 
F 7 * 77r (Q 2 ) and fit the CLEO data over the probed momentum range, denoted by {Ql xp }- In 
our recent analysis [§| the evolution ^(x; Q 2 ) = f/(ply — * Q 2 )^n(x] p| Y ) was performed for 
every individual point Q 2 xp , with the aim to return to the normalization scale p| Y and to 
extract the DA parameters (02, 04) at this reference scale for the sake of comparison with the 



previous SY results [6j. Stated differently, for every measurement, {Q 2 xp , ^ 7 ' 7 * 77r (Qexp)}> its 
own factorization and renormalization scheme was used so that the NLO radiative corrections 
were taken into account in a complete way. The accuracy of this procedure is still limited 
mainly owing to the uncertainties of the twist-4 scale parameter [9|, k ■ 5 2 , where the factor 
k expresses the deviation of the twist-4 DAs from their asymptotic shapes. (Another source 
of uncertainty originates from the mixingof the NLO approximations for the leading twist 
with the twist-4 contribution at LO, see [9J].) 

To summarize, the focal points of our procedure of the CLEO data processing are (i) 
a s (Q 2 ) is the exact solution of the 2-loop RG equation with the threshold M q = m q taken 
at the quark mass m q , rather than adopting the approximate popular expression in [l5| 
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that was used in the SY analysis. This is particularly important in the low-energy region 
Q 2 ~ 1 GeV 2 , where the difference between these two couplings reaches about 20%. (ii) All 
logarithms ln(Q 2 //i 2 ) appearing in the coefficient function are absorbed into the evolution 
of the pion DA, performed separately at each experimental point Q 2 xp . (hi) The value of 
the parameter 5 2 has been re- estimated in g to read 5 2 (lGeV 2 ) = 0.19 ± 0.02 GeV 2 . The 
present study differs from the SY approach in all these points and extends our recent analysis 
[9(| with respect to points (i) and (iii) yielding to significant improvements of the results. 

It turns out that the effect of varying the shapes of the twist-4 DAs exerts a quite strong 
influence which entails k to deviate from 1, i.e., from the asymptotic form. Note that next- 
to-leading-order corrections in the conformal spin |l6| for the twist-4 DAs cancel out exactly 
in the final expression for FX^" and therefore this deviation is due to more delicate effects. 
In the absence of reliable information on higher twists, one may assume that this uncertainty 
is of the same order as that for the leading-twist case. Therefore we set k — 1 ± 0.1. As a 
result, the final (rather conservative) accuracy estimate for the twist-4 scale parameter can 
be expressed in terms of k ■ <5 2 (lGeV 2 ) = 0.19 ± 0.04 GeV 2 , a value close to 0.20 GeV 2 used 
in Q. 

3. Confrontation with the CLEO data 

Pion DA vs the experimental data. To produce the complete 2a- and la-contours, 
corresponding to these uncertainties, we need to unite a number of regions, resulting from 
the processing of the CLEO data at different values of the scale parameter k ■ S 2 within this 
admissible range. This is discussed in technical detail in Here, we only want to emphasize 
that our contours are more stretched relative to the SY ones. The obtained results for the 
asymptotic DA (♦), the BMS model{K) 0, the CZ DA (■), the SY best-fit point (•10, a 
recent transverse lattice result (T) (l7l |. and two instanton-based models, viz., (*) [18[ and 
(♦) (using in this latter case m q = 325 MeV, n = 2, and A = 1 GeV) [l9j, are compiled in 
Tabled for the maximal, middle, and minimal twist-4 scale parameter. 



Table 1: Models/fits for different values of k ■ 5 2 (see text). 



k-5 2 0.23 GeV 2 


0.19 GeV 2 


0.15 GeV 2 


Models/fits (0,2,0,4,) 2 X 2 


(«2,«4) u 2 X 2 


(02,04) 2 X 2 


best-fit (+0.28,-0.29) 0.47 

• (+0.19,-0.14) 1.0 
X (+0.14,-0.09) 1.7 

♦ (-0.003, +0.00) 5.9 
■ (+0.40,-0.004) 4.0 
T (+0.06, +0.01) 3.8 

★ (+0.03, +0.005) 4.7 

♦ (+0.06,-0.01) 3.6 


(+0.22,-0.22) 0.47 
(+0.19,-0.14) 0.56 
(+0.14,-0.09) 0.89 
(-0.003, +0.00) 3.9 
(+0.40,-0.004) 5.2 
(+0.06, +0.01) 2.3 
(+0.03, +0.005) 2.9 
(+0.06,-0.01) 2.1 


(+0.16,-0.16) 0.47 
(+0.19,-0.14) 0.57 
(+0.14,-0.09) 0.52 
(-0.003, +0.00) 2.3 
(+0.40,-0.004) 7.0 
(+0.06, +0.01) 1.2 
(+0.03, +0.005) 1.6 
(+0.06,-0.01) 1.1 



We turn now to the important topic of whether or not the set of CLEO data is consistent 
with the nonlocal QCD sum-rule results for <p n . We present in Fig. ^ the results of the 
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Figure 1: Analysis of the CLEO data on F 7r7 * 7 (Q 2 ) in terms of error regions around the best-fit point 
(+) (broken line: la; solid line: 2a; dashed-dotted line: 3a) in the (112,0.4) plane contrasted with various 
theoretical models explained in the text. The slanted shaded rectangle represents the constraints on (02, a a) 
posed by the nonlocal QCD sum rules for the value A 2 = 0.4 GeV 2 . All constraints are evaluated at 
Uly = 5.76 GeV 2 after NLO ERBL evolution. 



data analysis for the twist-4 scale parameter k ■ 5 2 varied in the interval [0.15 < k ■ 5 2 < 
0.23] GeV 2 that includes both kinds of the discussed uncertainties of twist-4. We have 
already established in |?} that a two-parameter model <p n (x; 02, 04) factually enables us to 
fit all the moment constraints that result from nonlocal QCD sum rules (see jl^ for more 
details). It should be stressed, however, that the restriction on the Gegenbauer harmonics 
of order 2 is not just a plausible hypothesis but the direct result of the nonlocal QCD sum- 
rule approach for the pion DA. The next higher Gegenbauer harmonics up to the calculated 
order 10 turn out to be too small Q and are therefore neglected. The only parameter 
entering the nonlocal QCD sum rules is the correlation scale A? in the QCD vacuum, known 
from nonperturbative calculations and lattice simulations pH 22| . A whole "bunch" of 
admissible pion DAs resulting from the nonlocal QCD sum-rule analysis associated with 
A 2 = 0.4 GeV 2 at ^ 2 « 1 GeV 2 was determined [7|, with the optimal one given analytically 

by ^ MS {x) = < s (x) [l + a° pt • Cl /2 (2x - 1) + a° pt • C 4 3/2 (2x - 1)1 , where < s (x) = 6x(l - x) 

and a°i pt = 0.188, a° pt = -0.13 are the corresponding Gegenbauer coefficients. From Fig. 
E]we observe that the nonlocal QCD sum-rule constraints, encoded in the slanted shaded 
rectangle, are in rather good overall agreement with the CLEO data at the lex-level. This 
agreement could eventually be further improved by adopting still smaller values of A 2 , say, 
0.3 GeV 2 , which however are not supported by the QCD sum-rule method and also lattice 
calculations 2^1 . On the other hand, as it was demonstrated in the agreement between 



QCD sum rules and CLEO data fails for larger values of A 2 , e. g., 0.5 GeV 2 . 

Reliability of the main conclusions. The main qualitative conclusion of the presented 
analysis is that the "bunch" of pion DAs, derived in [?| from nonlocal QCD sum rules, agrees 
rather well with the CLEO data at the la-level, while both the CZ model and the asymptotic 
DA are ruled out at least at the 2<r-level. However, the value of the twist-4 contribution 
turns out to be a subtle point of the CLEO data processing. Having this in mind, let us 
inspect the stability of the main conclusions under the scope of the uncertainties associated 
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with this contribution. 

We have included the twist-4 parameter 5 2 in conjunction with the vacuum quark vir- 
tually 0, 5 2 ~ Aq/2. Now let us ignore this relation and assume that the total twist-4 
uncertainty is put by hand to an extreme uncertainty of, say, 30%, shifting the value of 
k ■ S 2 , at the low limit of the uncertainty, to k ■ S 2 = 0.13 GeV 2 . Would this change our 
conclusions dramatically? The result of this exercise is presented in Fig. Efa): One observes 
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Figure 2: Analysis of the CLEO data: (a) Assuming a twist-4 uncertainty of 30% or equivalently at 
S 2 = 0.13 GeV 2 . (b) Excluding the lowest 6 experimental points up to Q 2 — 3 GeV 2 . The designations here 
are the same as in Fig.^and the reference scale is /Xg Y = 5.76 GeV 2 . 

from Fig. EJa) that the asymptotic DA (♦) is still outside the 3a error contour (dashed- 
dotted line) and that the CZ point (■) is still far-away, while instanton-based models are 
just at the boundary of the 3a- (*) or 2a- (♦) ellipse. At the same time, the BMS model 
(K) moved practically to the center of the data region, whereas the Schmedding-Yakovlev 
best-fit point (•) ran outside the la-region. 

Another way to suppress the uncertainties of the twist-4 contribution is to repeat the 
processing of the CLEO data, excluding the low momentum transfer tail. At low Q 2 , the 
twist-4 contribution strongly affects the total form factor and, therefore, this exclusion can 
reduce the potential twist-4 uncertainties significantly. To study this effect in more detail, 
we removed in the data processing the lowest 6 experimental points (which possess very 
small errors) up to Q 2 xp = 3 GeV 2 reducing this way the relative influence of the twist-4 
contribution by factors of magnitude. Of course, the admissible cx-regions for the 02, 0,4 
parameters become much larger now due to this exclusion, as one sees from Fig. I21^b) in 
comparison with the LHS — a price one has to pay for the restricted way of data processing. 
Nevertheless, our main results and conclusions, discussed above, remain valid with the BMS 
model still inside the la ellipse and the asymptotic DA outside. 

The unknown high-order QCD radiative corrections provide another important source of 
systematic uncertainties. To estimate their size one should have at least the com plet e NNLO 
coefficient function of the process. A partial result, obtained quite recently in 23], gives a 
hint that the size of this contribution can be large. Therefore, the complete NNLO QCD 
calculation in the MS-scheme is a vital problem. In the absence of complete results, one 
can only roughly estimate the size of higher-order corrections by varying the reference scale 
fi 2 = n\ = fip = Q 2 , say, in the interval [Q 2 /2, 2Q 2 ]. The corresponding "shaking" of the 
form factor is taken into account in the systematic theoretical error that is demonstrated in 
Fig. |21 This uncertainty is rather large, of the order of la, and as a result, the set of the 
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Figure 3: Estimation of the influence of higher-order corrections by varying the reference scale in the range 
[Q 2 /2, 2Q 2 }. The designations are as in Fig.^ 

model predictions discussed above appears now inside or near the 2a contour (see Fig. OJ). 
Nevertheless our main conclusions remain valid. 

4. The inverse moment (x 1 ) 7r vs the CLEO data 

As already mentioned in the Introduction, in the present study we have processed the 
CLEO data in such a way as to obtain an experimental constraint on the value of the inverse 
moment (x _1 ) 7r ( / u 2 ) = f Q <p n (x; n 2 )x~ l dx that appears in different perturbative calculations 
of pion form factors. This is illustrated in Fig. |U^a), where the positions of the asymptotic 
DA, the CZ model, and the BMS one are also displayed. 

Fig. EJb) shows the theoretical estimate of the inverse moment obtained in the framework 
of nonlocal QCD sum rules. In fact, a "daughter sum rule" has been previously constructed 
directly for this quantity by integrating the RHS of the sum rule for ip n (x) with the weight 
x' 1 , (for details, see [tIHH)- Due to the smooth behavior of the nonlocal condensate at the 
end points x = 0,1, this integral is well defined eo ipso, supplying us with an independent 
QCD sum rule, with a rather good stability behavior of (a;~ 1 )| R (M 2 ), as one sees from 
this figure. Note that we have estimated (aT 1 )^/^ « 1 GeV 2 ) = 3.28 ± 0.31 at the value 
A 2 = 0.4 GeV 2 of the nonlocality parameter. It should be emphasized that this estimate is 
not related to the model pion DA, <y9 BMS (x; ci2, 04), constructed within the same framework. 
Nevertheless, the value obtained with the "daughter" QCD sum rule and those calculated 
using the "bunch" of pion DAs, mentioned above, (x _1 )^ MS (/xg) = 3.17 ±0.09 (jj, match 
each other. This fact provides further support for the self-consistency of the approach, as 
one appreciates by comparing the hatched strip with the BMS point (Jt) in Fig. Ufa). 

It is important to notice at this point that from the CLEO data one also obtains a 
constraint on the value of 02 + 04 = (x )^. xp (/i )/3 — 1 for the two Gegenbauer coefficients 
model. This constraint should be compared with the independent (from the theoretical 
model) estimate (x _1 )^ R (/io), as mentioned above. 

Let us discuss these results in more detail. In Fig.^a) we demonstrate the united regions, 
corresponding to the merger of the 2cr-contours (solid thick line) and the lcr-contours (thin 
dashed line), which have been obtained for values of the twist-4 scale parameter within 
the determined range (cf. Table [TJ. This resulting admissible region is strongly stretched 
along the (a 2 — a 4 ) axis, with the displayed models steered along (approximately) the same 
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Figure 4: (a) The result of the CLEO data processing for the quantity (:r _1 )° xp /3 — 1 at the scale 
Hq « 1 GeV 2 in comparison with the theoretical predictions from QCD sum rules, denoted SR. The thick 
solid-line contour corresponds to the union of 2cr-contours, while the thin dashed-line contour denotes the 
union of ler-contours. The light solid line with the hatched band indicates the mean value of (x- x )l R /3 - 1 
and its error bars in the second part of the Figure, (b) The inverse moment (x _1 )? R shown as a function 
of the Borel parameter M 2 from the nonlocal QCD sum rules at the same scale /Zg |7j ; the light solid line is 
the estimate for (x _1 )^ R ; finally, the dashed lines correspond to its error-bars. 



axis, demonstrating the poor accuracy for this combination of the DA parameters, while 
more restrictive constraints are obtained for (x^ 1 )^. One appreciates that the nonlocal 
QCD sum-rules result (x _1 )^ R , with its error bars, appears to be in good agreement with 
the constraints on (x _1 )° xp at the lcr-level, as one sees from the light solid line within 
the hatched band in Fig. Efa). In particular, the la-constraint obtained at the central 
value k ■ 5 2 = 0.19 GeV 2 exhibits the same good agreement with the corresponding sum- 
rule estimate because the theoretical uncertainty of the twist-4 scale parameter and of the 
radiative correction, already mentioned, affect mainly the ( ) constraint. The CLEO 

best-fit point (+) in Fig. Ufa) is near to zero in accordance with the previous data-processing 
results, presented in the first line of Table 1, 02 + cla — 0. Moreover, the estimate (x _1 )^ R is 
close to (x _1 ) EM /3 — 1 = 0.24 ± 0.16, obtained in the data analysis of the electromagnetic 
pion form factor within the framework of a different LCSR method in 2^, 26|. These three 
independent estimates are in good agreement to each other, giving robust support that the 
CLEO data processing, on one hand, and the theoretical calculations, on the other, are 
mutually consistent. Moreover, Dorokhov 1 recently obtained from the instanton-induced 
effective theory model <^ mod (a;) the estimate (ar~ 1 )^ nMd /3 - 1 « -0.09, which is close 
to the CLEO result. 

More importantly, the end-point contributions to the (x _1 )| R are suppressed, the range 
of suppression being controlled by the value of the parameter A 2 The larger this parameter, 
at fixed resolution scale M 2 > X 2 the stronger the suppression of the nonlocal-condensate 
contribution. Similarly, an excess of the value of (x^ 1 )^ over 3 (asymptotic DA) is also 
controlled by the value of A 2 , becoming smaller with increasing A 2 . Therefore, to match the 
value (x _1 )| R to the CLEO best-fit point, would ask to use larger values of X 2 than 0.4 GeV 2 . 
But this is in breach of the (a 2 , 04) error ellipses. A window of about 0.05 GeV 2 exists to 
vary A 2 : any smaller and one is at the odds with QCD sum rules and lattice calculations 
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2^1 ; any larger and the nonlocal QCD sum-rules rectangle can tumble out of the CLEO 
data region. 



5. Comparison with the E791 data 



Ver y re cently, an independent source of experimental data by the E791 Fermilab exper- 
iment [l2| has become available providing additional constraints on the shape of the pion 
DA. However, these data are affected by inherent uncertainties and their theoretical expla- 
nation by different groups 2^, 2^, 3(| is controversial. It is not our goal here to improve the 
theoretical framework for the calculation of diffractive di-jets diffraction. For our purposes 
it suffices to show basically two things: first, that our predictions for this process are not 
conflicting the E791 data and second, to show in comparison with other models for the pion 
DA that the BMS model has best agreement with these data, using for all considered models 
the same calculational framework. 
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Figure 5: Comparison of ip as (solid line), </? cz (dashed line), and the BMS "bunch" of pion DAs (strip, Q) 
with the E791 data 12]. The corresponding X 2 values are As: 12.56; CZ: 14.15; BMS: 10.96. 



To compare our model DA for the pion [7] with the E791 di-jet events [l2| and other pion 
DAs, we adopt the convolution approach developed in j^ij having also recourse to [U| . The 
results of the calculation are displayed in Fig. 03 making evident that, though the data from 
E791 are not that sensitive as to exclude other shapes for the pion DA (asymptotic and CZ 
model), also displayed for the sake of comparison, they are relatively in good agreement with 
our prediction. Especially, in the middle x region, where our DAs "bunch" has the largest 
uncertainties (see |7[), the predictions are not in conflict with the data. Note, however, that 
all theoretical predictions shown in Fig. El are not corrected for the detector acceptance. For 
a more precise comparison, this distortion must be taken into account. 



6. Conclusions 



Let us summarize our findings. They have been obtained by refining the CLEO data 
analysis, we initiated in jjj, in the following points. We corrected for the mass thresholds 
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in the running strong coupling and incorporated the variation of the twist-4 contribution 
more properly. In addition, the CLEO data were used to extract a direct constraint on the 
inverse moment (x _1 ) 7r (/ig) of the pion DA — at the core of form-factor calculations. This 
has relegated the CZ model and the asymptotic pion DAs beyond, at least, the 3cr-level 
(confidence level of 99.7%), with the SY best-fit point still belonging to the la deviation 
region (68%) in the parameter space of (02, 04), while providing compelling argument in favor 
of our model [jj, which is also within this error ellipse and remains there even assuming a 
potentially higher twist-uncertainty of the order of 30%. 

Both analyzed experimental data sets (CLEO ^ and Fermilab E791 01) converge to 
the conclusion that the pion DA is not everywhere a convex function, like the asymptotic 
one, but has instead two maxima with the end points (x = 0, 1) strongly suppressed — in 
contrast to the CZ DA. These two key dynamical features of the DA are both controlled 
by the QCD vacuum inverse correlation length Aq, whose value suggested by the CLEO 
data analysis here and in is approximately 0.4 GeV 2 in good compliance with the QCD 
sum-rule estimates and lattice computations 2^ . 
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